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Abstract 


The conventional pressure swing adsorption (PSA) is used as an alternative to cryogenic 
distillation. It has not yet proved to be competitive with cryogenic processes for large 
units. Moreover it yields only one component as pure product and the other being impure. 
Recently a new PSA, called as duplex PSA, has been reported. Recent studies indicate that 
it yields two products but sharp separation is not possible. In this work, a modified duplex 
PSA has been proposed to obtain two products with purities in excess of 99%. Simulation 
studies were carried out to assess the performance of the original and the proposed duplex 
PSA for the CO2-N2, CH4-H2 and N2-O2 systems. In the original duplex PSA, the 
maximum purity obtained of the heavy components (CO2 and CH4) was 90% and that of 
the light component (N2 and H2) was 95%. In the N2-O2 system, maximum purity of both 
the products obtained was 75%. In contrast, the purity of both the products improved 
dramatically to 99+% or more for all the systems studied in the proposed duplex PSA. 
Further, the productivity is higher than the conventional PSA. The proposed duplex PSA 
finds application in the capture and sequestration of CO2 from flue gases, upgradation of 
natural gas and as an alternative to cryogenic distillation of gas mixtures like air. 
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Chapter 1 


Introduction 

A highly energy intensive cryogenic distillation was the only option available for many 
years to fractionate air. The cryogenic distillation has now been supplemented by energy 
efficient and economical non-cryogenic technologies such as pressure swing adsorption 
(PSA) and membrane separation systems. PSA is the first commercial non-cryogenic 
process to produce pure nitrogen or oxygen from air. PSA has gained importance not 
only in air separation but also in the petrochemical industries. Process improvements and 
new generations of carbon molecular sieves and zeolites have resulted in increased 
efficiency of PSA systems since they were first introduced for commercial applications 
about 25 years ago. Around 5000 PSA systems have been installed worldwide over the 
last two decades, indicating the tremendous growth of the technology (Keith, 1997). 

Most of the conventional PSA processes developed so far are based on the Skarstrom 
cycle. These require a number of beds and large pressure ratio (Ph/Pl) leading to a large 
capital and operating costs. As an alternate to the Skarstrom PSA cycles and its variants, 
a new type of PSA cycle has been proposed by Leavitt (1992) and Hirose (1994). These 
are referred to as dual-reflux or duplex PSA cycle. The feed in a duplex PSA is 
introduced at some intermediate point along the length of the bed unlike Skarstrom type 
PSA where the feed is introduced at either ends of the bed. There is an inherent recycle of 
both the rich stream as well as the lean stream. 

In the conventional PSA, there is a narrow moss transfer zone which transverses along 
the length of the bed. The bed behind and ahead of the mass transfer zone is in 
equilibrium. This implies that at any time only a fraction of the bed is utilized for mass 
transfer. Whereas in the duplex PSA, mass transfer zone is extended throughout the bed 
length resulting in 5-10 times large productivity than that of the conventional PSA. 

One key feature distinguishing the duplex PSA from the conventional PSA is the 
concentration profile in the bed. In the duplex PSA, the concentration profiles in the bed 
are not destroyed and reformed during each cycle unlike in the Skarstrom type PSA. 
They are only shifted marginally depending upon the pressure ratio that is employed. 
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To assess the performance of the duplex PSA, simulation studies have been carried out 
with systems like carbon dioxide-nitrogen (CO2-N2) with nitrogen as non-adsorbing, 
methane-hydrogen (CH4-H2) with hydrogen as non-adsorbing and nitrogen-Oxygen (N2- 
O2) where both the components are adsorbing. The Original duplex PSA was not able to 
give clean separation. We proposed two modifications to the Original duplex cycle based 
on the simulation studies. They are referred to as MV-I duplex and MV-II duplex PSA 
cycle. 

The MV-I duplex, where the time of extract withdrawl has been modified, shows 
marginal improvement in the product purity over that of the Original duplex. Though it 
was a marginal improvement, it gave an insight to further modify the duplex cycle. In the 
MV-II duplex, where the step of extract withdrawl has been modified, the purity of both 
the products improved dramatically to 99 +% with 99 +% recoveries. 

The parametric studies such as effect of feed inlet position, raffinate recycle ratio, gas 
flow rate and desorption pressure have been done to assess the comparative performance 
of the Original duplex, MV-I duplex and MV-II duplex PSA. The raffinate recycle ratio 
is defined as the ratio of amount of gas recycled to the amount of raffinate product drawn. 
The MV-II duplex PSA seems to have potential to fractionate any binary gas mixture into 
clean products. 

Based on the analogy with the duplex PSA, another process with temperature swing, 
duplex TSA has been proposed and the possibility of separation has been explored with 
CO2-N2 system (N2 is non-adsorbing). 
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1.1 Structure of the thesis 


Chapter 2 presents origin and development of the duplex PSA, literature survey on 
duplex PSA and simulation studies on PSA. 

Chapter 3 explains the process description of the duplex PSA and duplex TSA 

Chapter 4 described mathematical modeling and method of simulation of duplex PSA 
and duplex TSA. 

Chapter 5 discusses the results of the simulation studies performed with duplex PSA and 
duplex TSA. 

Chapter 6 presents the conclusions of the study and also presents the scope for future 
work. 


3 



Chapter 2 


Origin and Development of duplex PSA 

The introduction of PSA processes is commonly attributed to Skarstrom and Guerin de 
Montgareuil and Domine in 1957-1958 (Ruthven, 1993). The Air Liquide process, 
developd by Montgareuil and Domine, utilized a vacuum swing, whereas the Esso 
process, pioneered by Skarstrom, used a low pressure purge to clean the adsorbent bed 
following the blowdown step. The first three applications of PSA namely air separation, 
air-drying and hydrogen purification were in fact foreseen and demonstrated by 
Skarstrom. These remain the most important practical applications for this technology, 
although newer processes such as carbon dioxide recovery and natural gas purification 
are gaining increased acceptance. 

The Skarstrom cycle in its basic form utilizes two packed adsorbent beds and four steps, 
(1) pressurization, (2) high pressure adsorption, (3) blowdown and (4) low pressure 
purge. Both beds undergo these four steps in a cycle. The Skarstrom cycle was only 
viable to get less selectively absorbable coihponent as pure product and generally not 
applicable for production of more selectively adsorbed component. To get the more 
selectively adsorbed component as product, Wilson (1982) developed a process that is 
superficially the inverse of the Skarstrom cycle. The Wilson cycle has four steps viz. (1) 
pressurization (2) purge at high pressure, (3) blowdown and (4) low pressure feed. 

The PSA processes which are based on Skarstrom cycle are limited by a single product 
with low recovery. Sirkar (1977) and Knaebel (1989) modified the Skarstrom cycle by 
providing reflux of both the products to yield high purities of both the products at 
reasonably high recoveries. 

A new type of PSA cycle has been proposed by Leavitt (1992) and Hirose (1994). These 
are referred to as duplex or dual-reflux PSA cycle. These cycles are characterized by feed 
at some intermediate point along the length of the bed unlike Skarstrom cycle and its 
variants which employ feed from end of the bed. Hirose et. al. (1992) has explained the 
principle behind the development of the dual-reflux PSA cycle and is explained below. 
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2.1 Duplex PSA: Origin and development 

Generally, the objective of PSA processes is either to remove or to concentrate one or 
many components from a binary or multicomponent feed gas mixture. The conventional 
PSA scheme for the purpose of removal (or stripping) is illustrated in Figure 2.1, which 
shows the bed conditions during a half cycle (of Skarstrom cycle). The feed gas with 
mole fraction Yf enters from the top of the high pressure bed and separated into a lean gas 
that leaves from bottom. A part of the lean gas is recycled to the low pressure bed as 
stripping recycle. The enriched gas leaves the low pressure bed from top. The mole 
fraction of heavy component in the lean gas, Ft can be lowered to a value nearly equal to 
0 when the operating conditions are selected properly. However, the mole fraction of the 
enriched gas. Ye caimot be increased at the same time up to a value higher than the 
pressure ratio (Ph/Pl)- This restriction is a thermodynamic limitation because the upper 
limit is controlled by the adsorption isotherm, a thermodynamic limit. Thus, a 
conventional PSA process, utilizing stripping recycle as the mode of separation is 
suitable to produce only light component as pure. 


Feed Enriched gas 



Lean gas 

5r«^'6 (R,Yr) 

Fig. 2.1 PSA process with stripping recycle 
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The conventional PSA scheme for the purpose of concentration (or enriching) is 
illustrated in Figure 2.2, which shows the bed conditions during a half cycle (of Wilson 
cycle). The feed gas with mole fraction Yf enters from the bottom of the low pressure 
bed. A part of the enriched gas from top of the low pressure bed is recycled to the high 
pressure bed as enriching recycle after compression. The lean gas leaves the high 
pressure bed from the top. The mole fraction of heavy component in the enriched gas. Ye 
can be increased to a value nearly equal to 1 when the operating conditions are selected 
properly. However, the mole fraction of the lean gas, Yr can not be lowered at the same 
time below a factor of the pressure ratio (Pl/Ph) because of the thermodynamic 
limitation. Thus, a conventional PSA process, utilizing enriching recycle as the mode of 
separation is suitable to produce only heavy component as pure. 


Enriched gas 
(E,Ye) 


Stripping limitation Lean gas 
YR>Yp(P,ypH) (R,Yr) 



Feed 

(F,Yp) 


Fig. 2.2 PSA process with enriching recycle 


The above two limitations can be overcome by reprocessing enriched gas leaving the 
stripping PSA in the enriching PSA, and lean gas leaving the enriching PSA in the 
stripping PSA. This concept leads to a new PSA process as shown in Figure 2.3 and 
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called as duplex or dual-reflux PSA process. It is a combination of the two conventional 
PSA processes represented by Figure 2.1 and 2.2. The feed gas is supplied at an 
intermediate position of the bed. By analogy to distillation, the feed inlet position divides 
each column into enriching and stripping sections, both at the same pressure. The length 
of the enriching and stripping section .depends on the feed position. For an appropriate 
combination of operating variables, this process can give two types of high purity 
products- both enriched gas and letin gas. 


Feed 

(F.Yk) 


Enriched gas 
(E, Ye) 



Lean gas 
(R.Yr) 


Fig. 2.3 Duplex PSA process with intermediate feed inlet position 
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2.2 Literature review on duplex PSA 

Though duplex PSA has been proposed in early 1990s, not much work is found in open 
literature. Hirose et. al. (1994) studied the removal and enrichment of carbon dioxide 
from the air-COa mixtures. Zeolite 1 3X was used as the adsorbent and parametric studies 
such as effect of feed inlet position, reflux ratio, gas flow rates, and pressure ratio on the 
enriched product concentration or recovery as well as lean product concentration were 
experimentally investigated (Hirose et. al. 1995). They found that the optimum feed inlet 
position was at about middle of the bed. The first numerical analysis on duplex PSA was 
carried out by Hirose et. al.(1996) for simultaneous removal and concentration of carbon 
dioxide dilute gas from air. They developed a simple isothermal model with negligible 
axial dispersion and pressure drops through the bed to investigate the effects of various 
combinations of operating parameters. Ritter et. al. (2002) has carried out experimental 
studies for separation and concentration of a dilute feed of ethane in nitrogen using 
WestVaco BAX- 1500 Activated carbon. The enrichment of useful trace components in 
air, such as Ne, Kr, Xe, etc. was studied by Ritter et. al. (2003) using only enriching 
reflux and parallel equalization. Equilibrium theory analysis of duplex PSA for separation 
of a binary mixture has been done by Ebner and Ritter (2004). 

Molecular gate process described in Keller patent (1982) uses pistons to produce cyclic 
gas flow and pressure variations at the two end of the beds, while feed is admitted at an 
intermediate point. This process has not been commercialized because they are difficult 
to scale up to commercialized size (Leavitt 1992). Farooq et. al. (2002) studied the 
performance of piston PSA process with feed at the end of the bed unlike Molecular gate 
process where feed was at intermediate point of the bed. 

2.3 Numerical simulation studies on PSA 

Most of the early simulation studies were carried out using the equilibrium model (Hill 
and Wong (1981); Knaebel and Hill (1985)). The equilibrium model provides valuable 
information about the dynamic behavior of the adsorption processes. Raghavan and 
Ruthven (1987) presented numerical simulation of PSA in the recovery of trace 
adsorbable from an inert carrier using the linear equilibrium and linear rate expressions. 
They also assumed that during pressure changes (blowdown and pressurization) there 
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was no mass transfer between the fluid and adsorbent. Raghavan and Ruthven discussed 
numerical simulation for a simple two bed PSA process in which effects of kinetics and 
changes in flow rate due to adsorption are significant. Alpay and Scott (1992) used the 
linear driving force (LDF) model to describe the adsorption and desorption kinetics in 
spherical particles. Chahbani and Tondeur (2000) discussed the mass transfer kinetics in 
PSA. They showed that the choice of the pore diffusion model plays a key role for 
obtaining reliable simulations. They compared their results with other models like LDF 
and equilibrium model. Kvamsdal and Hertzberg (1995) showed the effect of mass 
transfer during the blowdown step. According to their results, the frozen solid assumption 
is valid only in certain cases, and taking into account the mass transfer during the 
blowdown step gives a better overall model performance in the studied cases. Meunier 
and Huberson et al. (1992) studied a general non-isothermal model describing a single- 
component gas adsorption during pressurization with two mass intraparticle diffusions 
(macropore-micropore) and axial flow in the column. They found that the adsorption 
kinetics are not always a secondary effect and that they play an important role in the 
column dynamics as well as do the adsorption equilibrium properties. Lu and Rodrigues 
(1993) simulated a three- step one-column pressure swing adsorption process in 
isothermal and adiabatic cases for separation of a binary mixture. The bed dynaunics at 
the cyclic steady state have been studied in the three individual steps (pressurization, feed 
and blowdown) by looking at the axial profiles and histories of the process parameters at 
four positions inside the bed. Modeling PSA leads to a nonlinear dynamic partial 
differential equation system, which is difficult to solve. The parameters in PSA models 
affecting the behavior of the process are highly coupled to each other. The model 
equations were solved by orthogonal collocation and by finite difference methods by 
Raghavan and Ruthven (1985). Mate and Knaebel 1988) employed method of 
characteristics and found the effects of incomplete purge in PSA, Sun and Meunier 
(1991) introduced a solution adaptive gridding technique (SAG) combined with a four 
point quadratic upstream differencing scheme to satisfactorily resolve very sharp 
concentration and temperature variations occurring in the ease of small dispersing effects. 
Wen and Chou (1995) simulated a PSA process for air-5 A zeolite system using the 
Movii^ Finite Element Method (MFEM) using the property of moving nodes in MFEM 
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with small weighting factors. The use of small weighting factors corresponds to the 
strong action of the penalty function. They discussed the effects of the weighting factors 
and boundary conditions on the implementation of the MFEM. Sun et al. (1996) 
suggested an improved finite-difference method for simulation of diffusion-limited PSA 
process. They concluded that this method of solution adaptive gridding was well suited 
for describing the sharp transitions as it exhibits much less oscillations eind false 
diffusion. Morbidelli et al. (1997) used the orthogonal collocations on moving finite 
elements for simulating the PSA process. Hirose et. al. (2002) used a simple isothermal 
model for dual reflux PSA with negligible axial dispersion and pressure drop through the 
bed and validated it with experimental results. 
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Chapter 3 


Process Description: Duplex PSA and TSA 


This chapter presents the description of the Original duplex PSA, the MV-I duplex PSA, 
the MV-II duplex PSA and the duplex TSA. These are cyclic processes and are caixied 
out in two identical beds. The beds are cylindrical vessel packed with adsorbent particles. 
The adsorbent preferentially adsorbs or retains one of components (heavy) over the other 
(light). The amount adsorbed increases with pressure and decreases with temperature. 
The process gets reversed (i.e. amount gets desorbed) on decreasing the pressure or 
increasing the temperature. So, by swinging the pressure (pressure swing) or temperature 
(temperature swing) inside the beds and arranging them in a sequence of steps, a gas 
mixture can be fractionated into the clean products. 

The process details of the Original duplex, the MV-I duplex, the MV-II duplex PSA and 
the duplex TSA is given below. 

3.1 Original duplex PSA 

Original duplex PSA process consists of four steps in a cycle. Figure 3.1 shows the steps 
involved. 

Step 1: In this step, the feed is introduced into bed-1 and both the end products are 
drawn simultaneously (raffinate from bottom of the bed- land extract from top of the bed- 
2). Table below shows the status of the beds. 



Bed-1 

Bed-2 

Pressure 

' Ph 

Pl 

Status 

Feed 

Purge 


Feed is introduced at an intermediate point along the length of bed-1. The heavy 
component gets adsorbed inside the bed. A part of the outgoing stream, rich in light 
component is drawn from the bottom as raffinate product. The remaining part is recycled 
to the bottom of bed-2 as raffinate recycle. This will purge off the heavy component from 
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bed-2. A part of the outgoing stream from the top of bed-2, rich in heavy component is 
drawn as extract product. The remaining part is recycled to the top of bed-1 as extract 
recycle. The extract recycle is at low pressure, Pl so it has to be compressed to high 
pressure, Ph before recycling it into bed-1. 

Step 2 : This is pressure equalization and pressurization step. Table below shows the 
status of the beds. 



Bed-1 

Bed-2 

Pressure 

Ph ^Pl 

Pl ^Ph 

Status 

Blowdown 

Pressurization 


First, bottom end of both the beds are connected through a valve till pressure in the beds 
equalizes. Then, pressure in bed-1 is reduced to Pl and the gas coming out is used to 
pressurize bed-2 to Ph with the aid of a vacuum pump or compressor or both. It could be 
done from either ends of the beds viz. rich end (top of the beds enriched with heavy 
component) or lean end (bottom of the beds enriched with light component). In the next 
step, feed is switched to bed-2. 

These two steps complete the first half cycle. The other half cycle is the mirror image of 
the first half cycle. Step-3 is the mirror image of step-1 and step-4 is the mirror image of 
step-2. 

3.2 Modified version-I duplex PSA (MV-I duplex PSA) 

MV-I duplex PSA process has the same steps as that of the Original duplex PSA with only 
modification in the time of extract withdrawl. In Original duplex PSA, a part of the 
outgoing stream from the bed undergoing purge is collected as extract product 
continuously throughout the step time whereas in MV-I duplex PSA, first part of the 
outgoing stream from the bed undergoing purge is collected as extract product and the 
remaining amount is recycled to the top of the bed undergoing feed step as extract 
recycle. The purity of both the end products in the MV-1 duplex PSA gets enhanced over 
the Original duplex PSA. 
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Pressure history in bed: Pressure history in the beds in a cycle is the same for both 
the Original duplex and the MV-I duplex PSA. Figure 3.2 shows the pressure history in 
the beds in a cycle. In step 1 and step 3, pressure in the beds (undergoing feed and purge) 
does not change with time while in step 2 and step 4, pressure in the beds (undergoing 
equalization and pressurization) changes with time between Ph and Pl. 



Fig3.2 Pressure history in the bed in a cycle in original duplex and MV-I duplex PSA 


3.3 Modified version-II duplex PSA (MV-II duplex PSA) 

MV-II duplex PSA process consists of six steps in a cycle. Figure 3.3 shows the steps 
involved. In the Original duplex PSA both the end products are drawn simultaneously in 
a step whereas in the MV-II duplex PSA, the end products are dravm in the two different 
steps. 

Step 1: In this step, the feed is introduced into bed-1 and' only raffinate product is 
drawn. Table below shows the status of the beds. 



Bed-1 

Bed-2 

Pressure 

Ph 

Pl 

Status 

Feed 

Purge 
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t eed is introduced at an intermediate point along the length of bed-1. The heavy 
component gets adsorbed inside the bed. A part of the outgoing stream, rich in light 
component is drawn from the bottom as raffinate product. The remaining is recycled to 
the bottom ot bed-2 as raffinate recycle. This will purge off the heavy component from 
bed-2. I he total outgoing stream from the top of bed-2 rich in heavy component is 
recycled to the top of bed- 1 as extract recycle. The extract recycle is at low pressure, Pl 
so it has to be compressed to high pressure, Ph before recycling it into bed-1. 


Step 2: I'his is pressure equalization and pressurization step. Table below' shows the 
status of the beds. 



Bed-1 

Bed-2 

Pressure 

Ph ‘=t>Pi 

Pl ^Ph 

Status 

Intermediate-blowdown 

Pressurization 


hirst, bottom end of both the beds are connected through a valve till pressure in the beds 
equalizes. Then, pressure in bed-1 is reduced to Pi and the amount coming out is used to 
pressurize bed-2 to Ph with the aid of vacuum pump or compressor or both. But this 
amount is not sufficient to pressurize bed-2 to Ph as Pt < Pi so, extra amount is required. 
The extra amount is taken from the raffinate product drawn from bottom of bed- 1 during 
the step 2. By reducing the pressure in bed-1, the light component gets swept away from 
the bottom and the upper half portion becomes enriched with the heavy component. 


Step 3: This is extract withdrawl step. Table below shows the status of the beds. 



Bed-1 

Bed-2 

Pressure 

1 

Pi‘=^Pl 

Ph 

Status 

Final-blowdown 

Frozen 


Pressure in bed-1 is lowered from top end with the aid of a vacuum pump. The outgoing 
stream, rich in heavy component is collected as extract product. Final-blowdown 
pressure, Pl in bed-1 is governed by the intermediate-blowdown pressure. Pi and the 
amount of extract to be withdrawn. During this step, bed-2 is assumed to be frozen at Ph 
as time step is very small. It is assumed that during frozen step, there is no change in gas 
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phase or solid phase concentration profile with time in bed as frozen step time is small. 
These three steps complete the first half cycle. The other half cycle is the mirror image of 
the first half cycle. Step 4 is the mirror image of step 1, step 5 is that of step 2 and step 6 
is that of step 3. 

Pressure history in bed: Figure 3.4 shows the Pressure history in the beds in a cycle. 
In step 1 and step 4, pressure in the beds (undergoing feed and purge) does not change 
with time while in step 2 and step 5, pressure in the beds (undergoing equalization and 
pressurization) changes with time. In step 4 and step 6, one bed remains at constant 
pressure (frozen) while pressure in other bed (undergoing extract withdrawl or final- 
blowdown) changes with time. 



Fig.3.4 Pressure history in the bed in a cycle in MV-II duplex PSA 
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3.4 Duplex TSA 

Based on the analogy with the duplex PSA, another process with temperature swing, 
duplex TSA has been proposed and the possibility of separation has been explored. The 
duplex TSA process consists of two steps in a cycle. Figure 3.5 shows the steps. . 


E, Ye 





Step: 1 


E,Ye 



Step: 2 



Fig.3.5 Sequence of steps in duplex TSA in a cycle 


Step 1: In this step, the feed is introduced into bed-1 and both the end products are 
drawn simultaneously (raffinate from bottom of the bed- land extract from top of the bed- 
2). The operating pressure in both the beds are same. Table below shows the status of the 
beds. 



Bed-1 

Bed-2 

Status 

Cooling 

Heating 


Feed is introduced at an intermediate point along the length of the bed-1. Bed-1 is 
undergoing cooling so the heavy component gets adsorbed inside the bed. A part of the 
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outgoing stream, rich in light component is drawn from the bottom as raffinate product. 
The remaining part is recycled to the bottom of bed-2 as raffinate recycle. As bed-2 is 
undergoing heating, heavy component will get desorbed and recycle stream will purge off 
the desorbed heavy component. A part of the outgoing stream, rich with heavy 
component from the top of bed-2 is -drawn as extract product. The remaining part is 
recycled to the top of bed- 1 as extract recycle. 

Step-2 is the mirror image of step-1. In this step, feed is switched to bed-2. These two 
steps complete a cycle. 

Temperature history in bed: Temperature profile inside the beds has been assumed to be 
linear (minimum at the top and maximum at the bottom). During any step, one bed will 
undergo cooling and other bed will undergo heating and in next step cooling and heating 
will get reversed. Temperature shift between the two ends of the beds has also been taken 
as linear. This is the ideal case, and has been taken to study the possible mechanism of 
separation which will help in developing a real duplex TSA cycle. Figure 3.6 and 3.7 
shows the temperature shift in the beds undergoing cooling and heating respectively. 



Fig.3.6 Temperature shift along the length of the bed undergoing cooling 
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Chapter 4 


Modeling and Simulation 

This chapter presents the mathematical model and the method of simulation of the duplex 
PSA and the duplex TSA cycle. 

4.1 Mathematical model 

A gaseous mixture of N components is flowing through a fixed-bed of length L packed 
with adsorbent particles of identical size as shown in Figure 4.1. The flow is considered 
to be axial dispersed plug flow i.e. without gradients in the radial direction. The solid 
(adsorbent) and the fluid phases are considered to be separated by an interface. Mass and 
heat transfer between these two phases takes place through this interface. The governing 
partial differential equations are derived by making mass and heat balance in a 
differential element, Az of the fixed-bed adsorber as shown below. 



Fig 4. 1 Differential element of a section of the adsorber 


Mass balance for the bulk fluid phase 

The mass balance of component i over the differential element, Az for the gas phase is 
given as 

Input “ Output I- Accumulation 
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sAA-D^^ 


dz 


+ sA(vfC,)[ YsA[vj.c)\^^^^+6A^^z + A{\-s)^^z 


dc, . , ^ 5 ( 7 , 








J \~ 


Dispersion Convection 

where. 


"Y : Y- 

Dispersion Convection 


+A: dt 

J V. 


Accumulation 


Cj - gas phase molar concentration of component i, (mol/m^) 

Vf - interstitial fluid velocity, (m/s) 

A - cross sectional area of the adsorber , (m^) 
e - bed voidage 

Qi - solid phase molar concentration of component i , (mol/m^) 

t - time, (s) 

- molar rate of interfacial transfer of component i, (mol/m^.s) 
D, - axial dispersion coefficient, (m^/s) 


(4.1) 


Equation (4.1) can be rearranged to. 


eA 


/' dc. ^ 


V 


dz 


J 


^ dc '' 
-D,-^ + v,c, 


r+Ar 


V 


dz 


J 


+ Az + .4 (l - 6-)— Az = 0 

dt ^ ^ dt 


(4.2) 


On dividing Equation (4.2) by sAAz, we get 


dc, d ( \ 

dt dz^ ^ ’ 


dz 


dc.^ 

A.— 

- dz ) 


+(!z£)^=o 

s dt 


(4.3) 


Assuming the gas obeys ideal gas law, c, can be expressed as 

P 


c, = 


' RT. 


(4.4) 
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where, 


c, - total gas phase molar concentration (mol/m^) 

P - pressure {atm) 

R - universal gas constant {atm-m3/gmol-K) 

Tg - gas phase temperature (K) 

The gas phase concentration of component i is given by 

where, yj is the mole fraction of component i in the gas phase. 


On substituting the Equation (4.5) in Equation (4.3), we get 


^iC,y,) + -^{v,C,y,)-- 
ot oz oz 


D,- 


aC,y,y (1-£)8?, 


dz 


dt 


The pressure effect on axial dispersion coefficient can be expressed as 
D, 

i. p 

where Dio is the axial dispersion coefficient at reference pressure (1 atm). 


The interstitial fluid velocity, is calculated by the Blake-Kozeny equation. 




= -K 


dy dP 
I50ju{\-ef dz 

dP 


where. 


K 


dz 


p 


150//(l-ff) 


d - particle diameter, (m) 


On substituting Equation (4.5), (4.7) and (4.8) in Equation (4.6), we get 


(4.5) 


(4.6) 


(4.7) 


(4.8) 
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a Ur I az 


, 0 -^)% 0 

a^ierj p azUrJ z a 


On simplifying Equation (4.9), we get 


^'y' I I ^ ^ 

dt P dz^ P dz- P dz dz p-ydz) T; dz^ 7 ; dz dz PT^ dz dz T; ^ 5z 

I az' dz dz ? I a^ J T.dzdz p I a/ J 2 ; I, a/ J p s dt 


for i -1 to N-1 


(4.10) 


On summing up Equation (4.10) over all components X ~ ^ rearranging, we get 


ap D,.„ a-p ifapV. p d% dpdr^ 2P ( dT^ 

dt~ p dz^ plazj T. dz^ T. dz dz r„M & 


,a^P rSPV PdPdT^ pfdp 


dz^ [dz T dz dz T\dt) ^ s dt 


V 


(4.11) 


« dP 

Further simplification of ~Qt^~dt 
Equation (4. 1 0) can be rearranged as, 

— +^— ^ 
dt [ P^ J dz^ dz dz P dz dz 


dz^ dz dz dz dz ^ P^ J P V a? 


^ 4 k-^ 


yJdP' 


P dz^ dz- dz dz P/ ( dz 


P 1, ar J T^[dt 


RTA^-^)dq, 

P e dt 


(4.12) 


If 


K±^ 
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Equation (4.12) simplifies further to. 


ay. _ D„ dy, y, d‘T, 2 si; , 2y, ( dT, 


dt P 5z“ dz dz T^" y dz J J 

dz~ dz dz T dz dz p[dz J p[dt J t\ dt j P s dt 


On summing up Equation (4.13) over all components ' and rearranging, we get 


dt J Tg 52^ I ydz) dz dz \ T^[ dt 


e M dt 


(4.14) 


By putting Equation (4.14) into Equation (4.13), further reduces to 



2 dT^ 

Lf;3y,eP^RT,{i-e)f 

N 

y.Y. 

dq, 

.fiLl 

^dz^ 

dz dz 

j dz dz P e V 

/=i 

dt 

dt y 


for i =7 to N-1 


(4.15) 


For the isothermal process, 

dT d^T 

Tg = constant ■=> — - = 0 ; — y- = 0 , so Equation (4.14) and (4. 1 5) will reduce to 

dz dz 


dt dz^ [dz ) ^ s dt 


(4.16) 


dy< _D,Jd^y>'\ ^ dP ^ RT, {^-s) f ^dq, dq, 
dt P \ dz^ \ dz dz P el ' /=i dt dt 


for i -1 to N-1 (4.17) 


25 


Mass balance for the solid phase 

The rate of interfacial mass transfer between the solid and the fluid phase is represented 
by the Glueckauf s linear driving force (LDF) model. 


dt 




for i -1 to N 


(4.18) 


where, 

kf - linear driving force(LDF) mass transfer coefficient, (s'‘) 

q] - solid phase molar concentration in equilibrium with c. , {mol/m^) 

Heat balance for the bulk fluid phase 

Considering the process is adiabatic, the heat balance over the differential element, Aa 
shown in the Figure 4.1 for the gas phase is given by 

Input = Output + Accumulation 


sA\-K 


T" 

Conduction 


dz 




= 6A{-K 


e^ax 


dz 


Convection 




Conduction 


z+Az 




Y 

Convection 


V ^ ^ 

Interfacial transfer 




(4.19) 


V 

Accumulation 


One assumption has been made that the velocity, of the fluid is constant in the small 
differential element Az and will vary from the element to element along the length of the 
bed. 

On dividing Equation (4.19) by Az and rearranging, we get 


dt~ pf^dz^ ^ dz ^ 


(4.20) 


26 



The velocity is calculated by the Blake-Kozeny Equation (4.8). On putting Equation (4.8) 
into Equation (4.20), we get 


dt 

where. 


= a- 


d% 


+ K 


dz dz 


n^a. 








^c,ax ~ effective axial thermal conductivity of gas, (W / mK) 
p - density of gas, {kg hrP) 

- specific heat of gas, ( J / kgK ) 

K 

(X . — ^ effective axial thermal diffusivity of gas, {m^/s) 

Pfp. 

hp - particle heat transfer coefficient, (W/m^K) 

- density of solid, {kg/n?) 

C^, - specific heat of solid, ( J / kgK ) 

- solid phase temperature, (K) 


Heat balance for the solid phase 

Considering the process is adiabatic, the heat balance over the differential element, Az 
shown in the Figure 4.1 for the solid phase is given by 

Input = Output + Accumulation 
X(-A//),;^(l-^)^Az + V,(7; -2 ;)^Az = |;[^(1-£)aC,,X]Az 

/s=i ut Ol 

'' ^ ^ ^ '''' ^ 

Heat of adsorption Interfacial transfer Accumulation 


(4.21) 


(4.22) 
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where, 


a 




; specific surface area of the adsorbent particle, (m ' ) 


AH - heat of adsorption, (J/mol) 


On dividing Equation (4.22) by A([-s)Az and rearranging, we get 



1 




/ = ! 


SQ, I 

ACp,v(l-£) 




( 4 . 23 ) 


The duplex PSA cycle has been simulated at isothermal condition whereas the duplex 
TSA cycle has been simulated at non-isothermal condition. 

The governing partial differential equations (PDE’s) derived above for duplex PSA and 
duplex TSA are being summarized below. 


Governing PDE’s for duplex PSA 


dt 


K 


dz^ I dz 


^ dt 


(4.16) 


A- 

dt P 




dP (1 - g) r ^ dq^ dq^ 
dz dz P e dt dt 


for i=ltoN-l (4.17) 


dt 


kfqj-q,) 


for i =i to N 


(4.18) 
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Governing PDE’s for duplex TSA 


dt 


Du 




Z: 


dz 


-RZ. 


(l-g) ^ dq. 

e it dt 


Z„ dz^ 


+ K 


p^^fdPy PdPdZ^ 


dz- 


\dz j 


Z^ dz dz 


( dZ. ^ 


+ - 


Z.. 


\ dt J 


( 4 . 14 ) 


dt P 


r ^2 


aV, 2 dy^ dZ^ 


dz~ Z^ dz dz 


+ K 


dy. dP ^ RZg (l ^) ( i 
dz dz P g ' /=] dt dt 

for i =1 to N-1 


(4.15) 


dt 


for i -1 to N 


( 4 . 18 ) 


ar a^r_ bp dz„ h a , ^ 

— L = a—f + K—-^ tL±-(T -t\ 


dz dz g p^C 


PS 


(4.21) 


ar 1 




+- 




a "a p.c,(x-c) 




(4.23) 


4.2 Initial and boundary conditions 

The partial differential Equations [(4.14), (4.15), (4.16), (4.17) and (4.21)] are first order 
with respect to time and second order with re.spect to length. It requires initial condition 
and two boundary conditions whereas Equations [(4.18) and (4.23)] are first order with 
respect to time and requires only initial condition to solve it numerically. The initial and 
boundary conditions for different steps of duplex PSA and duplex TSA are given below. 
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4.2.1 Original duplex and MV-I duplex PSA 

The Original duplex PSA and the MV-I duplex PSA cycle employs four steps viz. 
blowdown, purge, pressurization and feed step. The initial and boundary conditions for 
each step are given below. 


> Blowdown step 
I.C’s: att=0 

yi (z, 0) = jv/, ; P (z, 0) = Pn 
qi(z,0) = q*i(y,P,T^ 



> Purge step 

I.C’s: at t =0 

yi (z. 0) =yi (z, ts) ; P (z, 0) = P (z.ts) 

qi(z,0) = qi(z,tB) 



& r=0 


K 



= 0 
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> Pressurization step 

l.C's: at t =0 

y, (z, 0) =yi (z, tp) ; P (z, 0) ^ P (z,tp) 
q, (z, 0) = qi (z, tp) 



> Feed step 

l.C’s: att=0 

yi (z, 0) = ;;; (z, tpp) ; P (z, 0) = P (z, tpp) 
qi(z,0) = qi(z,tpp) 



z,r and z,* represents the bed position just before and after the feed inlet position 
respectively. Boundary conditions at these two grids are given below. Expression for 
composition of the mixed stream, y^„ has been given in Appendix-A. 
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5 ;^ 


dz 


= 0 


dt 


= 0 


at z,. 


dz 




dz 


IL 

K 


4.2.2 MV-II duplex PSA 

The MV-II duplex PSA cycle employs six steps viz. intermediate-blowdown, final- 
blowdown, purge, pressurization, frozen and feed step 

There are two blowdown steps viz. intermediate-blowdown and final-blowdown step. 
The boundary conditions for the intermediate-blowdown and the final-blowdown step are 
same as that of the blowdown step in the Original duplex PSA. Other steps like purge, 
pressurization and feed steps have the same boundary conditions as that of the Original 
duplex. 

4.2.3 Duplex TSA 

The duplex TSA cycle employs two steps viz. heating and cooling step. The initial and 
boundary conditions for each step are given below. 

> Heating step 
l.C's: att=0 

yi(z,0) =yuni; Tg(z, 0) =Tg,ini 

Ts (z, 0) = Ts,i„i ; qi(z, 0) = g*i(yi. P, Tg) 
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ar„ 


dz 


a 


(T. 


«■»> 


) 

^ z=0 


dT„ 


dz 


= 0 




> Cooling (feed) step 

1. C 's: at t =0 

yi(z,0) =-yi(z,th); 

Ts (z, 0) = Ts (z, ti) ; 


Tg(z, 0) = Tg(z, ti) 
qi (z, 0) = qi (z, ti) 



dz 

dz 

dT 


Z*bO 


A. 


(y,„ 


rsO 


= -^ 
K 


dz 


a 


(T ■ -T\ ^ 

^ g,"’ g I -=0 ^ 


dz 

dt 

dT^ 

dz 


= 0 


z=L 


= 0 




zf and zf" represents the bed position just before and after the feed inlet position 
respectively. Boundary conditions at these two grids are given below. Expression for 
composition and temperature of the mixed stream, y-„ and T^ j„ is given in Appendix-A. 


at zf 


dz 

dt 

ar 


= 0 


= 0 


at Zp'^ 


dz 


Z=Zf 


Z=Zf 


dz 


0 


dz 

dZ 


dz 


Z^Zf 


.Z/ 

K 

Vf . I 

C(T , -t\ ) 


Z=Zp 


a 
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4.3 Parameters estimation 

The estimation of parameters used in the mathematical model of the duplex PSA and the 
duplex TSA is described below. 


Adsorption equilibrium isotherm 

The equilibrium solid phase concentration has been represented . by the extended- 
Langmuir isotherm 



b,c, 


/=! 


for i = 1 to N 


where, 

j - saturation constant of component i , (mol/m^) 
bf - Langmuir constant of component i, (mol/m^)"' 
temperature dependency of qsi and bi (Cho et al., 1995) 

exp(A:ir^,) 

expC^/g) 


(4.24) 


b,=b^^exp{k^T) 

where 

ki, k 2 - temperature parameters, (K) 
Tg - gas phase temperature, (K) 


Interfacial mass transfer coefficient 

The interfacial mass transfer coefficient has been estimated by the Glueckauf model. 




r: 


(4.25) 


<lif 


where, 

ki - linear driving force mass transfer coefficient, {s'^) 
Sp - particle porosity 

- effective diffusivity, (mVj) 


34 



Rp - particle radius, (m) 

Cif - feed concentration, {mol/m^) 

qjj - solid phase molar concentration in equilibrium with feed concentration 
of component I (mol/m^) 


The validity of the LDF model however is limited to 


D t 

e C 

K 


>0.1 but it can also provide 


an adequate description of the diffusion behaviour for even 


— ^ < 0.1 , provided that the 


LDF constant is adjusted as a function of 


R} 


(Alpay and Scott, 1992). 


where tc is half cycle time, (s) for a cyclic process. 


Axial dispersion coefficient 

The axial dispersion coefficient has been estimated from the correlation given by Wen 
and Fan (1975) 


where 


O.SdpVf 

D,,,=0.75A. + ^ 


'M 


1+9.5 


\^p^f J 


is molecular diffusivity, (m^/s) 


(4.26) 


Particle heat transfer coefficient 

The particle heat transfer coefficient has been estimated from the correlation given by 
Wakao and Kaguei, (1982) 




fiC 




pg 


K. 


g J 


2.876 0.3023 
Re, Re;^^ 


(4.27) 


where, 

hp = particle heat transfer coefficient, (Pf / m^ K) 
Gg = Pgff ; mass velocity of fluid, (Kg/m^s) 
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Re 


p 


fJ- 


; particle Reynolds number 


Effective axial thermal conductivity of gas 

The effective axial thermal conductivity for the gas flowing through a packed bed has 
been estimated from the correlation given by Wakao and Kaguei, (1982) 

^=|i+0.5(Pr)(Re,) (4,28) 

where, 

Ke, ax - effective axial thermal conductivity of gas 

Kg - thermal conductivity of gas 

Ke° - effective thermal conductivity for quiescent bed 

Pr - Prandtl number 

Rcp - particle Reynolds number 


Kf for a bed of spherical particle is given by 








m = 0.280 - 0.757 log,,, s - 0.057 log 
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where, 

Ks - thermal conductivity of solid 
£ - bed voidage 


(4.29) 


Energy equation 

Energy required for the compression or evacuation of a gas stream in a open system is 


given by (Smith and Van Ness, 1996) 


W (isentropic) = 


ynRP 


y-l 


/ D ^ 


Y-l 


y 


-1 


(4.30) 
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where, 


y = Cp/ Cy = 1 .67 (monoatomic gas) 

1 .40 (diatomic gas) 

1 .30 (polyatomic gas) 

n - moles of gas to be compressesed or evacuated from pressure P i to P 2 
/ = 1 .4 has been taken in energy calculation. 


Valve equation 

The valve equation is based on a simple pipe flow and equivalent length corresponding to 
the pressure drop through a valve 


P, 

On rearranging the above equation, we get 




(4.31) 


where 


v,=c(APr 


C = 


^ Dgc 


P-fl-eP 
f - fanning friction factor 
- density of gas, (kg/m"^) 

Ig - equivalent length of pipe, 
D - diameter of the pipe, (m) 


In the blowdown and the pressurization step, velocity at the opening of the bed should be 
such that the bed reaches the desired pressure in the given step time. For this, the above 
equation has been modified as 

Vf=kC{LPf" 

The k is found by the regula-falsi method. 


37 



4.4 Method of solution 

The governing partial differential equations [(4.19), (4.20) and (4.21)] for the duplex PSA 
and [(4.22), (4.23), (4.24), (4.25) and (4.26)] for the duplex TSA have been 
simultaneously solved along with their initial and boundary conditions. To solve them 
numerically, the PDE’s were discretized in space domain using central difference scheme 
of the finite difference method (FDM). Thus the PDE’s were converted into a set of 
ordinary differential equations. The discretized form of the PDE’s for both the duplex 
PSA and the duplex TSA is given in Appendix-A. Since the equation is stiff in nature, it 
has been solved using an ODE solver, namely DDASPG subroutine in the IMSL library 
of FORTRAN-PowerStation which is based on the Gear’s method. 

The simulation code has been written in FORTRAN sind executed on a Pentium-4 PC 
with 1GB RAM. CPU time was found to depend mainly upon the number of grids used in 
discretization of the PDE’s. CPU time taken for one cycle with 53, 105, 157 and 205 
grids are 5, 17, 38 and 77 s respectively. The cyclic steady state was reached between 150 
to 200 cycles depending on the operating parameters. 

Flow charts for simulation of the Original duplex PSA and the MV-I duplex PSA and the 
MV-II duplex PSA cycle is given below in Figure (4.2) and (4.3). 
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START 
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Specs: F, Y,-, E, /?«, T,, P^, P,., L, d, t,. tp, tpp, t,, icyde"'^^ 
Data: Langmuir isotherm parameters (b„ qd 


Initialise: P = P„, y,=yt^ .j 


icyde = 1 
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Set inlet velocity; 


sAP,t,, 


Solve PDE’s to get V/, q,, P profile in the bed 


Calculate final amount in the bed 


Gear 's 
Algorithm 




Calculate from mass balance 

(i.e. output = inlet + initial amount - final amount) 


Read tn, Ph 



Find new v/ by 
regula falsi method 


\Ph-P\ 

within ti>n 


Calculate final amount in the bed 



Calculate from mass balance 

(i. e. input = final amount - initial amount) 




























Read ty, F and F/ 


Set final value of qt, P from pressurization as initial condition of feed 


Calculate initial amount in the bed 


Newton-Cotes 

formula 


Set reflux inlet velocity; v 
Set feed velocity; v. .. 


icvcle + 1 


Solve PDE’s to gety/, qi, P profile in the bed 


Algorithm 


X 

Calculate final amount in the bed 
1 

^ 

4 

Newton-Cotes 

formula 



Calculate balance 

(Le. output = reflux inlet + F + initial amount -final amount) 


icvcle = icycle 


Set final value of>;„ q,, P 
from feed as initial 
condition of blowdown 


























Final-blo wdown Intermediate-blo wdown 


START 



Specs: F, Y/,-, E, Tp Pff, Pi, L, d, tm, tm.-, tp, tpp, tp-, 
Data: Langmuir isotherm parameters (bp q^i) 


Initialise: P = yi=yi q, = qi*(Pn, v, 


icycle = 1 


Read tm, Pi 


Calculate initial amount in the bed 


Newton-Cotes 

formula 


Assume Vf 


Solve PDE’s to getjj, qi, P profile in the bed 


Gear 's 
Algorithm 


Find new vy by 
regula falsi method 


P, -P\<i 
within /«/ 


Calculate final amount in the bed 


Newton-Cotes 

formula 


Calculate from mass balance 

(i.e. output = initial amount - final amount) 


Read tsF 


Set final value ofy* q,, P from intermediate blowdown as initial condition of final blowdown 


Calculate initial amount in the bed 


Newton-Cotes 

formula 


Assume Pl 


Assume Vf 
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Feed Pressurization 
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Calculate n'^'p , from mass balance 
(i.e. output - reflux inlet + + initial amount— final amount) 











Chapter 5 


Results and Discussion 

This chapter presents the results of the simulation studies of the Original duplex PSA, the 
MV-1 duplex PSA, the MV-II duplex PSA and the duplex TSA. 

The following set of systems has been used to assess the comparative performance of the 
three duplex PSA cycle. 

> Systems with one component as non-adsorbing. 

• Carbon dioxide and nitrogen over Zeolite 13X (nitrogen is non-adsorbing). 

• Methane and hydrogen over Activated carbon (hydrogen is non-adsorbing). 

> System where both the components are adsorbing. 

• Nitrogen and oxygen over Zeolite 5A 

The possibility of separation in the duplex TSA cycle has been presented based on the 
analogy with the duplex PSA. 

The product purity, recovery, productivity and energy requirement has been used as 
measures to assess the process performance. These are defined as; 


Product purity (%) = 


moles of desired component in product 
Total moles of product 


X 100 


Recovery (%) = 


moles of desired component in product 
moles of desired component in feed 


X 100 


Productivity = 


moles of de sired component in product 
(weight of adsorbent) (time) 


Energy requirement has been calculated for: 

• Compression of a- flowing stream from the low pressure bed to the high pressure 
bed 

• Pressurizing a bed from low to high pressure and 

• Evacuating a bed from high to low pressure. 
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5.1 Fractionation of carbon dioxide-nitrogen (CO 2 -N 2 ) mixture 

Carbon dioxide is one of the major green house gases causing global warming. The flue 
gases discharged by power plants and industries contribute to more than 70% of total CO 2 
emitted into the atmosphere. Therefore it is important to capture CO 2 from flue gases to 
maintain eco-balance of the earth. Currently, absorption of CO 2 using various solvents is 
being practiced commercially. The energy required for the solvent regeneration is 
significantly high. Adsorptive and membrane separations are the potential alternatives to 
the existing absorption processes. We have assessed the performance of the duplex PSA 
for recovery of CO 2 from the flue gas. In the present work, we have taken flue gas as a 
binary mixture of CO 2 and N 2 . As N 2 adsorbed inside the bed is very small compared to 
that of the CO 2 , we have considered N 2 as non-adsorbing component. The parameters 
used in simulation are listed in Table 5.1. 

Table 5.1 Parameters used in simulation of fractionation of carbon dioxide-nitrogen mixture 


1 . Molar feed composition 

20% CO 2 , 80% N 2 

2. Feed position 

middle of the bed 

3. Adsorbent 

Zeolite 13X 

4. Particle diameter, cm 

0.0707 

5. Bed voidage 

0.4 

6. Bulk density, kg/m^ 

590 

7. Bed length, cm 

100 

8. Bed diameter, cm 

2.5 

9. Langmuir constant for CO 2 (b), m^/mol 

0.4405 

10. Saturation constant for CO 2 (qs), mol/m^ 

3411.0 

1 1 . LDF constant for CO 2 , s'^ 

0.808 

12. Operating temperature, K 

298.15 

1 3 . Cycle time, s (Original and MV-I duplex) 

50 

Feed (tp), s 

20 

Blowdown (ts), s 

5 

Purge (tp), 5 

20 

Pressurization (tpa), s 

5 
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14. Cycle time, 5 (MV-II duplex) : 54 

Feed (tp), 5 - . 20 

Intermediate-blowdown (tsi), 5 : 5 

Final -blowdown (tBp), ^ ; 2 

Purge (tp), ^ : 20 

Pressurization (tpp), 5 : 5 

Frozen (tpR), 5 : 2 


The equilibrium data for carbon dioxide-nitrogen mixture over Zeolite 13X reported by 
Hirose et al. (1996) has been used. Figure 5.1 shows the adsorption isotherm of carbon 
dioxide on Zeolite 13X for different pressure at 298.15 K. 



Fig.5.1 Adsorption isotherm of carbon dioxide on Zeolite 13X at 298.15 K 


Simulation results of the Original duplex, MV-I duplex and MV-II duplex PSA is given 
below. 
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5,1.1 Original duplex PSA 

In Original duplex PSA, the equalization and pressurization step can be operated in two 
modes. It could be done either from the CO 2 rich end (mode-I) or the N 2 rich end (mode- 
II). Figure 5.2 shows two different modes of equalization and pressurization. 


CO 2 rich end 



N2 rich end 

Mode-I Mode-II 


Fig.5.2 Two modes of equalization and pressurization 

The performance of both the modes has been compared. Figure 5.3 and 5.4 shows the 
mass balance sheet for the mode-I and the mode-II respectively. 

The product purity is about the same in both the modes. But, the energy requirement in 
the mode-I is about 1.5 times more than that in the mode-II. This is because of the 
amount of gas to be compressed or evacuated is more in the mode-I than that in the 
mode-II. The amount of extract recycle to be compressed is about same in both the modes 
but the amount coming out from the bed undergoing blowdown in mode-I is about 7 
times more than the amount coming out in mode-II resulting in more energy requirement. 
In mode-II, most of the CO 2 gets locked inside and only a small amount of CO 2 comes 
out from the bed. We call it as ‘locking in blowdown’ and is explained below. 

Locking in blowdown: The amount adsorbed inside the bed is governed by its partial 
pressure in the gas phase. To understand the locking of CO 2 , we have examined and 
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compared the partial pressure and solid phase profile of CO2 in the bed undergoing 
blowdown for both the modes of equalization and pressurization. 

Figure 5.5 (a) shows the partial pressure and solid phase concentration profile of CO2 in 
the bed undergoing blowdown step from CO2 rich end (mode-I). The initial and final 
partial pressure profile of CO2 in the first half of the bed is almost same and there is 
significant change only in the second half. The same trend can be seen for the solid phase 
profile. The change in solid phase profile in second half of the bed is significant so a 
large amount of CO2 comes out from the bed. 

Whereas, on evacuating the bed from the N2 rich end, a small amount of CO2 comes out 
because of the locking effect. Figure 5.5 (b) shows the partial pressure and solid phase 
concentration profile of CO2 in the bed undergoing blowdown from N2 rich end (mode- 
II). It is clear from the solid phase concentration profile that the CO2 gets desorbed and 
then again adsorbed in the first half of the bed according to the partial pressure profile 
while, in the second half of the bed, there is no significant change in the profiles. The 
overall effect is that, a small amount of CO2 comes out from the bed or most of the CO2 
gets locked inside the bed. 
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It can be seen that, the product purity is about the same for both the modes but energy 
requirement is less in the mode-II. Hence, operating duplex PSA in mode-II is more 
effective than mode-I from energy consideration. So now onwards, we have considered 
the equalization and pressurization of the beds to be done from the N 2 rich end. 

Figures 5.6 (a, b, c, d, e, f, g, h) show the change in the gas and the solid phase 
concentration with time in the bed undergoing feed, blowdown, purge and pressurization 
(equalization and pressurization from N 2 rich end) for the mass balance sheet given in 
Figure 5.4. In all the profiles, dimensionless bed length, ‘0’ represents the end of the bed 
where the stream is entering. The direction of incoming and outgoing streams from the 
beds has been represented by arrow ( — ►) sign. 
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Feed: Figures 5,6 (a) and 5.6 (b) show the gas and the solid phase concentration profiles 
of the feed step. 1 here are two inlet streams. The extract recycle, rich in CO 2 is fed to the 
top end and the feed is introduced in the middle of the bed. The CO 2 gets adsorbed and 
the raffinate product, rich in N 2 is drawn from the bottom. At the feed position, sudden 
change in the bed profiles can be observed as feed mixes with a stream of different 
composition. The CO 2 present in the extract recycle gets adsorbed only in the initial 
portion of the bed and the latter portion up to the feed position seems to be unused. The 
bed profiles can be stretched by increasing the amount of the recycle but then product 
purity will go down. 



0.0 0.1 0.2 03 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Dimensionless bed length 


Fig,5.6 (a) Gas phase concentration profile in the bed undergoing feed step (Original duplex) 

3500 
3000 

I 

1 2500 

t'i 

0 

^ 2000 

1 

«u 

-g, 1000 
^ 500 

0 

0.0 0.1 0.2 03 0,4 0,5 0.6 0.7 0.8 0.9 1.0 

Dimensionless bed length 

Fig. 5.6 (b) Solid phase concentration profile in the bed undergoing feed step (Original duplex) 
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Blowdown: Figures 5,6 (c) and 5.6 (d) show the gas and the solid phase concentration 
profiles of the blowdown step. As the pressure reduces from Ph (1 atm) to Pl (0.3 atm), 
CO 2 gets desorbed from the solid and the gas composition of CO 2 increases. It can be 
seen that in the first half of the bed, CO 2 gets locked inside and most of the CO 2 comes 
out from the other half of the bed. The overall effect is that a very small amount of CO 2 
(0.64 mmol) comes out from the bed and the outgoing stream is enriched with N 2 . 



Fig.5.6 (c) Gas phase concentration profile in the bed undergoing blowdown step (Original duplex) 



Fig. 5,6 (d) Gas phase concentration profile in the bed undergoing blowdown step (Original duplex) 
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Purg€, I^igurts 5.6 (c) and 5.6 (f) show the gas and the solid phase concentration profiles 
ot the purge step. As the raffinate recycle stream rich in N2 enters the bed, the gas phase 
composition oi CO2 decreases. This results in the decrease in the partial pressure of CO2 
as total pressure is constant (0.3 atm). So, the CO2 gets desorbed and purged off from the 
bed. A part of the outgoing stream, rich, in CO2 is collected as extract product. 



0,0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Dimensionless bed length 


CO 2 enriched 
extract 


Fig.5.6 (e) Gas phase concentration profile in the bed undergoing purge step (Original duplex) 



‘ig.5.6 (f) Solid phase concentration profile in the bed undergoing purge step (Original duplex) 
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Pressurization: Figures 5.6 (g) and 5.6 (h) show the gas and the solid phase 
concentration profiles of the pressuiization step. The bed is pressurized from Pl (0.3 atm) 
to Ph (f utni) with the gas coming out (5.7 mmol) from the bed undergoing blowdown. 
But, the amount of gas lequired to pressurize the bed is 5.8 mmol. The extra amount (0.1 
mmol) is taken from the outgoing stream of the bed undergoing feed step. It can be seen 
from the bed piofiles that the composition of CO2 in the gas phase decreases but the 
change in the solid phase amount is marginal and only in the initial portion of the bed. 
This is because the inlet stream is rich in Nj and most of the CO2 present gets adsorbed 
on just entering the bed. 



Fig.5.6 (g) Gas phase concentration profile in the bed undergoing pressurization step (Original duplex) 



0.0 o.i 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Dimensionless bed length 

Fig.5.6 (h) Solid phaso eoncantration profile in the bed4irMlergoing pressurization step (Original duplex) 
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5.1.2 Parametric study of Original duplex PSA 

The effect of feed, raffinate recycle ratio and desorption pressure on the performance of 
Original duplex PSA has been presented below. 

Effect of raffinate recycle ratio, Rr 

Figure 5.7 shows the effect of raffinate recycle ratio on the product purity and energy 
requirement at two different values of feed rate. At a low value of Rr (0.1), product 
purity is small. This is because the amount of gas recycled is not sufficient enough to 
purge off the CO 2 from the bed undergoing purge. On further increasing the amount of 
gas recycled, product purity gets improved and reached an optimum at Rr value of 0.2. 
Further increase in Rr resulted in the decrease of product purity. The continuous decrease 
in purity is because of the non adsorbing nature of the N 2 . The bed (purge) in which the 
raffinate recycle enters is at low pressure 0.3 atm and at this pressure the void space 
inside the bed can hold only 2.6 mmol of the gas. So, most of the N 2 present in the 
recycle stream will come out from the bed resulting in poor performance of the Original 
duplex. It can be seen from Figure 5.7 that even at two different values of feed, the 
optimum purity is reached at Rr value of 0.2. 

On increasing the amount of raffinate recycle, the amount of gas circulating between the 
beds increases so the energy required to compress the extract recycle stream coming out 
from the low pressure bed also increases. 



•a 

(U 


o 

£ 


E 

53 

■£ 

CJ" 

(D 

>% 

00 


C 

m 


Fig. 5.7 Effect of raffinate recycle ratio; PH=latm, PL-0.3atm 
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Effect of feed, F 

Figure 5.8 shows the effect of feed on the product purity and the energy requirement. As 
the amount of feed increases, the amount of raffinate recycle stream (rich in N2) increases 
for a fixed value ol raffinate recycle ratio. As explained above, most of the N2 present in 
recycle stream will come out from the bed affecting the extract purity. 

As feed increases, the energy requirement per unit mole of feed decreases. This is 
because of the decrease in relative amount of the extract recycle (i.e. Re/F) to be 
compressed from Pl( 0.3 atm) to Ph (1 atm). 



10 30 50 70 90 110 130 150 170 

Feed (mmol/half cycle) 


Fig.5.8 Effect of feed; PH=latm, PL=0.3atm, Rr=0.2 

Effect of desorption pressure, Pl 

Figure 5.9 shows the effect of desorption pressure on the product purity and the energy 
requirement. At large value of desorption pressure (0.9 atm), purity of both the end 
products is poor. This is because at this pressure, CO2 is strongly adsorbed inside the bed 
and the raffinate recycle is not able to properly purge off the adsorbed CO2. Product 
purity increases on decreasing the desorption pressure and reaches an optimum value at 
about Pl value of 0.3 atm. On further decreasing the desorption pressure, concentration 
front of CO2 in the bed undergoing blowdown breaks through and more amount of CO2 
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starts coming out. As this stream is used for pressurization, the CO 2 present gets adsorbed 
around the end of the bed from where N 2 rich raffinate product will be drawn in the next 
step. Because of this, a slight drop in the product purity can be observed. 

As P[, decreases, the energy requirement increases because of increase in the pressure 
ratio (Ph/Pl)- 



Fig.5.9 Effect of desorption pressure; F=60mmol/half cycle, PH=latm, Rr=0.4 


From the above pstfametric study it can be concluded that the product purity in Original 
duplex is limited by the raffinate recycle gas (rich in non adsorbing N 2 ). To overcome the 
recycle limitation, we modified the time of extract withdrawl in Original duplex. The 
basis of the modification is explained below. 

5.1.3 MV-I duplex PSA 

In Original duplex, a part of the outgoing stream from the bed undergoing purge is 
collected as extract product (rich in CO 2 ) continuously throughout the step time (20 s). 
Figure 5.10 shows the variation in the mole fraction of CO 2 and corresponding 
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cumulative amount drawn with time in outgoing stream from the bed undergoing purge. 
It can be seen that the initial amount coming out is pure CO 2 and then, the CO 2 
composition starts decreasing and drops to 70% after 20 s. This is because of the N 2 
present in the gas (raffinate recycle) entering the bed starts coming out continuously with 
time. 



Fig.5.10 Variation in mole fraction of CO 2 and cumulative amount with time in outgoing 
stream from purge; F=100mmo 1/half cycle, PH=latm, PL=0.3atm, Rr=0.2 

So there is a scope to improve the product purity by modifying the time of extract 
withdrawal. The first part of the outgoing stream from the bed undergoing purge is 
collected as extract product and the remaining amount is recycled as extract recycle to the 
top of the bed undergoing feed step. By doing so, the purity of both the end products gets 
enhanced over the purities obtained with the Original duplex PSA. 

Figure 5.11 shows the mass balance sheet of the MV-I duplex Avith operating parameters 
same as that of the Original duplex (Figure 5.4). 
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It can be seen that composition of CO 2 in the rich gas (extract) improved from 77% to 
8 1 % and that of N 2 in lean gas (raffinate) from 94% to 95%. 

However, the trend in the gas and the solid phase concentration profile in the bed is the 
same as that of the Original duplex. 

The comparative performance of the MV-I duplex and the Original duplex PSA is given 
below. 

5.1.4 Comparison of MV-I duplex and Original duplex PSA 

Figure 5.12, 5.13 and 5.14 shows the comparative performance of the MV-I duplex and 
the Original duplex PSA with change in raffinate recycle ratio, feed and desorption 
pressure respectively. It can be seen that in all cases, performance of MV-I duplex 
improved over Original duplex. 

Figure 5.12 and 5.14 shows that at low value of raffinate recycle ratio (0.1) or large value 
of desorption pressure (>0.5) there is no improvement in product purity because CO 2 is 
not properly purged off. 

Effect of raffinate recycle ratio 



Fig.5.12 Effect of raffinate recycle ratio; F=60mmol/half cycle, PH=latm, PL=0.3atm 




Product purity (mole fraction of CO2) 







Though theie is a marginal improvement in product purity in MV-I duplex over Original 
duplex, it gives an insight to fuither modify the Original duplex cycle for clean separation 
ofC02-N2. 

5.1.5 MV-n duplex PSA 

From the parametric study of Original and MV-I duplex PSA, it can be concluded that the 
product purity is limited by the recycle stream circulating between the beds. The extract 
product cannot be enriched to 99+% CO 2 when the extract is drawn from the outgoing 
stream ol the bed undergoing purge step. To overcome this, the extract withdrawl step in 
MV-II duplex has been modified and it dramatically improved the product purities. 

It can be seen that during equalization and pressurization step, most of the Ni present in 
the bed undergoing blowdown gets swept away from the bottom and the upper half 
portion becomes enriched with pure CO 2 . On further evacuating the bed from top end, 
pure CO 2 as extract will come out. 

Figure 5.15 shows the mass balance sheet of MV-II duplex for 100 mmol/half cycle of 
feed. The adsorption pressure, Ph has been fixed as 1 atm and intermediate-blowdown 
pressure, Pj as 0.5 atm. The final-blowdovm pressure, Pl (0.34 atm) is governed by the 
intermediate-blowdown pressure and the amount of extract (20 mmol) to be drawn. 
During equalization and pressurization, the amount coming out from bed-1 tmdergoing 
intermediate-blowdown is not sufficient to pressurize bed-2 to Ph as Pl < Pi and an extra 
amount of 1 .4 mmol is required. The extra amount is taken from the raffinate product 
drawn from bottom of bed- 1 undergoing feed step. The purity of extract in MV-II duplex 
improved drastically from 77% to 99.5% with 99.7% recovery and that of raffinate from 
94% to 99.9% with 99.9% recovery. 

Figures 5.16 (a, b, c, d, e, f, g, h, i, j) show the change in the gas and the solid phase 
concentration with time in the bed undergoing feed, intermediate-blowdown, final- 
blowdown, purge and pressurization step for the mass balance sheet given in Figure 5.15. 
In all the profiles, dimensionless bed length, ‘0’ represents that end of the bed where the 
stream is entering. The direction of incoming and outgoing streams from the beds has 
been represented by arrow ( — ►)sign. 
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Fig.5.15 Mass balance sheet of MV-II duplex PSA 













Feed: Figures 5.16 (a) and 5.16 (b) show the gas and the solid phase concentration 
profiles of the feed step. There are two inlet streams. The total gas (rich in CO 2 ) coming 
out from the bed undergoing purge is fed to the top as extract recycle, and feed is 
introduced in the middle of the bed. The CO 2 gets adsorbed and the raffinate product 
(rich in N 2 ) is drawn from the bottom. At the feed position, sudden change in the bed 
profiles can be observed as the feed mixes with a stream of different composition. The 
bed profiles are stretched throughout the length as the amount of extract recycle (164 
mmol) is quite large compared to that in the Original duplex (42 mmol). 



Fig.5.16 (a) Gas phase concentration profile in the bed undergoing feed step (MV-II duplex) 



Fig.5.16 (b) Solid phase concentration profile in the bed undergoing feed step (MV-II duplex) 
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Intermediate-blowdown: Figures 5.16 (c) and 5.16 (d) show the gas and the solid phase 
concentration profiles of the intermediate-blowdown step. As the pressure reduces from 
Ph (1 atm) to P] (0.5 atm), the CO 2 gets desorbed from the solid and the gas composition 
of CO 2 increases. A little amount of CO 2 (0.25 mmol) comes out from the bed because of 
the locking effect. The N 2 present in first half of the bed gets swept away and becomes 
enriched with pure CO 2 . 



Fig.5.16 (c) Gas phase concentration profile in the bed undergoing intermediate-blowdown step (MV-II duple? 



Fig.5.16 (d) Solid phase concentration profile in the bed undergoing intermediate-blowdown step (MA^-Il duple 
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Final-blowdown: Figures 5.16 (e) and 5.16 (f) show the gas and the solid phase 
concentration profiles of the final-blowdown step. It can be seen from the bed profiles 
that initially, upper half portion of the bed is enriched with pure CO 2 . On evacuating the 
bed fiom top, the outgoing gas (extract product) is enriched with 99.48% CO 2 .To draw 
20 mmol extract product, bed has to be evacuated from P, (0.5 atm) to Pl (0.34 atm). 

CO 2 enriched 



Fig.5.16 (e) Gas phase concentration profile in the bed undergoing final-blowdown step (MV-II duplex) 



Fig.S.lb (f) Solid phase concentration profile in the bed undergoing final-blowdown step (MV-II duplex) 
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Purge: Figures 5.16 (g) and 5.16 (h) show the gas and the solid phase concentration 
profiles of the purge step. The bed undergoing purge is at final-desorption pressure (0.34 
atm). As the raffinate recycle stream rich in N 2 enters the bed, the gas phase composition 
of CO 2 decreases. This results in the decrease in the partial pressure of CO 2 as total 
pressure is constant. So, the CO 2 gets desorbed and purged off from the bed. 



Fig.5.16 (g) Gas phase concentration profile in the bed undergoing purge step (MV-II duplex) 



Fig.5.16 (h) Solid phase concentration profile in the bed undergoing purge step (MV-II duplex) 
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Pressurization: Figures 5.16 (i) and 5.16 (j) show the gas and the solid phase 
concentration profiles of the pressurization step. The bed is pressurized from Pi (0.34 
atm) to Ph (1 atm). It can be seen from the bed profiles that the composition of CO 2 in the 
gas phase decreases but the change in the solid phase amount is marginal and only in the 
initial portion of the bed. This is because the inlet stream is rich in 95% Ni which is non- 
adsorbing. 



Fig.5.16 (i) Gas phase concentration profile in the bed undergoing pressurization step (MV-II duplex) 



Fig.5.16 (j) Solid phase concentration profile in the bed undergoing pressurization step (MV-II duplex) 
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5.1.6 Parametric study of MV-II duplex PSA 

The effect of feed, raffinate recycle ratio and desorption pressure on the performance of 
MV-II duplex PSA has been studied. 

Effect of raffinate recycle ratio, Rr 

Figure 5.17 shows the effect of raffinate recycle ratio on the product purity and the 
energy requirement. At low value of Rr (0.2), product purity is small. On increasing the 
recycle ratio, product purity gets improved. There is no change in product purity between 
recycle ratio 0.4 to 1 .4 and thereafter that it starts decreasing. 

On increasing the raffinate recycle ratio, the amount of gas circulating between the beds 
increases so the energy requirement also increases. 



The reason for the purity trends can be understood by looking at the bed profiles of the 
intermediate-blowdown step at different recycle ratios. It gives the details of how much 
N 2 is left inside the bed at the end of this step. Since the bed after intermediate-blowdown 
step goes for withdrawl of CO 2 rich extract (final-blowdown step), the N 2 present inside 
the bed may come out with the CO 2 effecting its purity. Figures 5.18 (a), (b) and (c) 
show the gas phase concentration profiles in the bed undergoing intermediate-blowdown 
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for a recycle ratio of 0.2, 0.6 and 1.8 respectively. It can be seen that at a Rr value of 
0.2, the CO 2 in extract recycle entering from top of the bed is not sufficient and they get 
adsorbed only in the initial part. Most part of the bed at the end of feed step is filled with 
Ni. The N 2 present in the bed is not properly swept away even after evacuating the bed 
(during intermediate-blowdown). As. the bed after intermediate-blowdown goes for 
extract withdrawl, the N 2 left in the initial portion of the bed comes out with CO 2 
effecting its purity. Whereas, at raffinate recycle ratio of 0.6, the first half portion of the 
bed after intermediate-blowdown gets enriched with pure CO 2 and on evacuating the bed, 
extract (rich in CO 2 ) with 99.5% purity comes out. 

The drop in product purity at Rr value of 1.8 can be explained in the same way from the 
corresponding bed profile shown in Figure 5.18 (c). 

Effect of feed, F 

Figure 5.19 shows the effect of feed on the product purity and the energy requirement. 



20 60 100 140 180 ,220 260 

Feed (mtnol/half cycle) 


Fig.5.19 Effect of feed; PH=latm, Pi=0.5atm, Rr=0.6 

At a feed rate of 60 mmol/half cycle, the extract purity is 97% and that of raffinate is 
99%. On increasing the feed, the purity increases and reaches an optimum (extract purity 
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- 99.5% and raffinate purity - 99.7%) at a feed rate of 100 mmol/half cycle. On further 
increasing the feed rate, a slight drop in product purity is observed. The reason can be 
explained in same way as that of the raffinate recycle ratio effect because increasing the 
feed is equivalent to increasing the amount of recycle stream. 

The energy requirement increases with feed rate because of the decrease in Pl required to 
draw the corresponding amount of the extract (F*Xf). 

Effect of intermediate-desorption pressure. Pi 

Figure 5.20 shows the effect of intermediate-desorption pressure on the product purity 
and the energy requirement. At large value of Pi (> 0.4 atm), purity of both the end 
products are poor. This is because during equalization and pressurization step, desorption 
pressure (Pi) is not sufficient to sweep off the N 2 present in the bed undergoing 
intermediate-blowdown. The N 2 left inside the bed will come out when extract (rich in 
CO 2 ) will be drawn. Other reason is that at large value of Pi, the corresponding pressure 
(Pl) in the bed undergoing purge is also large and the raffinate recycle is not able to 
properly purge off the adsorbed CO 2 , as it is strongly adsorbed. On decreasing the 
intermediate-desorption pressure, product purity gets improved and at Pi value of 0.3 atm 
purity of both the end products reached 99.9%. 



Fig.5.20 Effect of intermediate-desorption pressure; F=60mml/half cycle, PH=la(m, Rr=0.4 
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As P| deci eases, the corresponding pressure (Pl) in the bed undergoing purge also 
decreases. As the outgoing stream from purge (extract recycle) has to be compressed 
hom P|, to I H (1 atm) before recycling it into the bed undergoing feed, the increase in 
energy requirement is obvious with decrease in Pi. 

5.1.7 Comparison of MV-II duplex and Original duplex PSA 

Figures 5.21, 5.22 and 5.23 show the comparative performance of MV-II duplex and 
Original duplex PSA with change in raffinate recycle ratio, feed and desorption pressure 
respectively. It can be seen that in all cases, the product purity in MV-II duplex improved 
drastically over the purities obtained with the Original duplex. 

Effect of raffinate recycle ratio 



Raffinate recycle ratio 

Fig.5.21 Effect of raffinate recycle ratio; F=100nimol/half cycle, PH=latm, 
Pi==0.5atm (MV-II), PL=0.3atm (Original) 
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Product purity (mole fraction of CO2) ^ Product purity (mole fraction of CO2) 


Effect of feed 



Feed (mmol/half cycle) 

Fig.5.22 Effect of feed; PH=latm, P,=0.5atm (MV-II), PL=0.3atm (Original) 




P] fox MV-II \ Pl for Original (atm) 


Fig.5.23 Effect of intermediate-desorption (MV-II) or desorption pressure (Original); 
F=60mmol/halfcycle, PH=latm, Rr=0.4 
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5.2 Fractionation of methane-hydrogen (CH 4 -H 2 ) mixture 

The separation of CH 4 -H 2 is important in petroleum and chemical process industries. We 
have assessed the performance of the duplex PSA for clean separation with CH4-H2 
mixture. The parameters used in simulation are listed in Table 5.2. 


Table 5.2 Parameters used in simulation of fractionation of methane-hydrogen mixture 


1 . Molar feed composition 

30% CH 4 , 70% Hi 

2. Adsorbent 

Activated carbon 

3. Particle diameter, cm 

0.0707 

4. Bed voidage 

0.4 

5. Bulk density, kg/rn^ 

480 

6 . Bed length, cm 

100 

7. Bed diameter, cm 

2.5 

8 . Langmuir constant for CH 4 (b), m^/mol 

3.91E-03 

9. Saturation constant for CH 4 (qs), mol/m 

5067.0 

10. LDF constant for CH 4 , s'^ 

3.01 

1 1 . Operating temperature, K 

298.15 

12. Cycle time, s (Original and MV-I duplex) 

50 

Feed (tp), s 

20 

Blowdown (tfi), s 

5 

Purge (tp), j 

20 

Pressurization (tpa), s 

5 

13. Cycle time, s (MV -II duplex) 

54 

Feed (tp), s 

20 

Intermediate-blowdown (tsi), s 

5 

Final-blowdown (tep), s 

2 

Purge (tp), 5 

20 

Pressurization (tpR), s 

5 

Frozen (tpa), s 

2 
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The equilibrium data for methane-hydrogen mixture over Activated carbon reported by 
Malek and Farooq (1996) has been used. The amount of hydrogen adsorbed is negligible 
over a wide range of pressure and it can be considered as non-adsorbing component 
Figure 5.24 shows the adsorption isotherm of methane on Activated carbon for different 
pressures at 298. 1 5 K. 



Fig.5.24 Adsorption isotherm of methane on Activated carbon at 298.15 K 

As composition of feed in the CH 4 -H 2 system (Xf= 0.3) is different from that of the CO 2 - 
N 2 system (Xp = 0.2), we first studied the effect of feed position on the product purity. In 
CO 2 -N 2 system, optimum feed position was middle of the bed (Hirose et. al. 1995). We 
observed that in CH 4 -H 2 system, optimum feed position has shifted towards the top 
(extract recycle end) of the bed. 

5.2.1 Effect of feed position 

Figure 5.25 shows the effect of feed position on the product purity. In the study, the total 
length of the bed has been divided into 53 equispaced grids. It can be seen that on 
introducing feed in the middle, the system performance is poor. The highest purity is 
obtained on introducing feed at 10* grid (10/53 grid) from extract recycle end of the 
bed. 
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Fig.5.25 Variation in product purity with different position of feed; F=60 mmol/half cycle, 

Ph=1 atm, Pl= 0.5 atm, Rr=0.5 



Fig.5.26 Final gas phase concentration profile of CFI4 in bed undergoing feed; F-60 mmol/half cycle, 

Ph= 1 atm, Pl= 0.5 atm, Rr=0.5 


The reason can be understood by analyzing the final gas phase concentration profile in 
the bed undergoing the feed step. Figure 5.26 shows the profile corresponding to different 
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feed location. On introducing the feed between 5*'’ to 15*’’ grids from extract recycle end, 
there is smooth transition in the profile throughout the bed length. Whereas, on 
introducing feed after 20 “^ grid from extract recycle end, feed breaks through the bed 
resulting in poor product purity. 

Here onwards we have taken 10* grid (optimum) from extract recycle end as the feed 
position (with total number of grids equal to 53 ). 

The comparative performance of the Original duplex, MV-I duplex and MV-II duplex 
PSA has been studied with variation in raffinate recycle ratio, feed and desorption 
pressure respectively. It has been observed that the trend in product purity with variation 
in simulation parameters is same as that of the CO 2 -N 2 system studied earlier. This is 
obvious as in both the systems, only one component gets adsorbed and other is non- 
adsorbing. 

5.2.2 Original duplex and MV-I duplex PSA 

Figure 5.27 shows the mass balance sheet of the Original duplex PSA. For a feed rate of 
60 mmol/half cycle, purity of extract (rich in CH 4 ) is 78.5% and that of raffinate (rich in 
H 2 ) is 90.3%. 

As H 2 is not getting adsorbed, it affects the extract purity in same way as that of N 2 (non- 
adsorbing) in CO 2 -N 2 system. 

Figure 5.28 shows the mass balance sheet of the MV-I duplex with operating pairameters 
same as that of the Original duplex (Figure 5.27). In MV-I duplex, the first part (I) of the 
outgoing stream from the bed undergoing purge is collected as extract product (rich in 
CH4) and the remaining amount (II) is recycled to the top of the bed undergoing feed as 
extract recycle. It can be seen that in MV-I duplex, the extract purity improved from 
78.5% to 87.9% and that of the raffinate from 90.3% to 94.4% over Original duplex. 
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System: CH4- Hz 

Adsorbent : Activated carbon Purity (%) 78.50 90.33 

Recovery (%) 78.5 90.3 

Productivity (mol/kg-h) 4.0 10.7 



duplex PSA 





































5.2.3 Comparison of MV-I duplex and Original duplex PSA 

Figures 5.29, 5.30 and 5.31 show the comparative performance of the MV-I duplex and 
the Original duplex PSA with change in raffinate recycle ratio, feed and desorption 
pressure respectively. It can be seen that in all cases, performance of MV-I duplex 
improved over Original duplex. The trend is similar to that of the CO 2 -N 2 system and can 
be explained in the same way. 

Figures 5.29 and 5.31 show that at low value of raffinate recycle ratio (0.2) or at a large 
value of desorption pressure (>0.5) there is no improvement in product purity because 
CH4 is not properly purged off from the bed undergoing purge step. 

Effect of raffinate recycle ratio, Rr 
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5.2.4 MV-II duplex PSA 

From the parametric study presented in section 5.2.3, it can be concluded that the clean 
separation is not feasible in Original and MV-I duplex PSA because of the non-adsorbing 
nature of H 2 . 

We studied the performance of MV-II duplex to get clean separation with CH 4 -H 2 
system. Figure 5.32 shows the mass balance sheet of MV-II duplex for 60 mmol/half 
cycle of feed. It can be seen that in MV-II duplex, the extract purity improved 
dramatically from 78.5% to 99.9% with 99.9% recovery and that of the raffinate from 
90.3% to 99.9% with 99.9% recovery over Original duplex. 

The gas and solid phase concentration profile of the bed undergoing intermediate- 
blowdown, final-blowdown, purge and pressurization step has not been shown as they 
exhibit almost the same trend as that of the CO 2 -N 2 system. There is change only in the 
feed profile as feed location is different. Figures 5.33 (a) and (b) show the gas and the 
solid phase concentration profiles in the bed undergoing feed step. 
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5.2.5 Parametric study of MV-II duplex PSA 

The effect of raffinate recycle ratio, feed and desorption pressure on the performance of 
MV-II duplex PSA has been studied. The trend in product purities and energy 
requirement is the same as that of the CO 2 -N 2 system. 

Effect of raffinate recycle ratio, Rr . 



Effect of feed, F 
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Fig.5.35 Effect of feed; PH=latm, Pi=0.5 atm, Rr-0.5 
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Effect of intermediate-desorption pressure, Pj 
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Fig.5.36 Effect of intermediate-desorption pressure; F=60 mmol/half cycle, Pn-latm, Rr-0.4 


5.2.6 Comparison of MV-II duplex and Original duplex PSA 

Figures 5.37, 5.38 and 5.39 show the comparative performance of MV-II duplex and 
Original duplex PSA with change in raffinate recycle ratio, feed and desorption pressure 
respectively. It can be seen that in all cases, the product purity in the MV-II duplex 
improved drastically over the purities obtained with the Original duplex. 
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Product purity (mole fraction of CH4) Product purity (mole fraction of CH4) 


Effect of raffinate recycle ratio 
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Raffinate recycle ratio 

Fig.5.37 Effect of raffinate recycle ratio; F=60 mmol/half cycle, PH=latm, 
Pi=0.6 atm (MV-II), Pl=0.4 atm (Original) 
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Fig.5.38 Effect of feed; PH=latm, P,=0.S atm (MV-II), Pl=0.4 atm (Original) 
Rr=0.5 (MV-II), Rr=0.4 (Original) 
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Product purity (mole fraction of CH4) 


Effect of desorption pressure 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Pi for MV-II; Pl for Original (atm) 

Fig.5.39 Effect of intermediate-desorption (MV-II) or desorption pressure (Original); 
F=60mmol/ha!f cycle, PH=latm, Rr=0.4 
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5.3 Fractionation of nitrogen-oxygen (N 2 -O 2 ) mixture 


The equimolar feed composition has been taken to study the fractionation of nitrogen- 
oxygen system. In this system, both the components get adsorbed unlike the case of CO 2 - 
N 2 and CH 4 -H 2 systems. We have explored the possibility of fractionation in the duplex 

PSA. The parameters used in simulation are listed in Table 5.3. 

Table 5.3 Parameters used in simulation of fractionation of nitrogen-oxygen mixture 

1 . Molar feed composition 

50% N 2 , 50% O 2 

2. Adsorbent 

Zeolite 5 A 

3. Particle diameter, cm 

0.0707 

4. Bed voidage 

0.4 

5. Bulk density, kg/rr^ 

670 

6. Bed length, cm 

100 

7. Bed diameter, cm 

2.5 

8. Langmuir constant for N 2 (b), m^ /mol 

2.813E-03 

9. Saturation constant for N 2 (qs), mol/m^ 

5260.0 

1 0. LDF constant for N 2 , s'^ 

19.7 

11 . Langmuir constant for O 2 (b), m^/mol 

8.935E-04 

12. Saturation constant for O 2 (qs), mol/rr^ 

5260.0 

13. LDF constant for O 2 , 

62.0 

14. Operating temperature, K 

298.15 

15. Cycle time, ^ (Original duplex) 

50 

Feed (tp), s 

20 

Blowdown (te), s 

5 

Purge (tp), 5 

20 

Pressurization (tpR), s 

5 

16. Cycle time, 5 (MV-II duplex) 

60 

Feed (tp), s 

20 

Intermediate-blowdown (tsi), 5 

5 

Final-blowdown (tsp), s 

5 

Purge (tp), 5 

: 20 


94 




Pressurization (Ipr), s 
Frozen (Ifr), 5 


5 

5 


The equilibrium data for nitrogen-oxygen mixture over Zeolite 5 A reported by Farooq et. 
al. (1996) has been used. Figure 5.40 (a) shows the adsorption equilibrium isotherm (total 
solid phase amount vs. gas phase composition of N 2 ) for different pressures at 298.15 K. 
Figure 5.40 (b) shows that the separation factor is 3.15 in favor of nitrogen. 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Gas phase compostion of N2 

Fig.5.40 (a) Adsorption isotherm of nitrogen-oxygen on Zeolite 5 A at 298.15 K 



Fig.5.40 (b) Solid phase composition vs. gas phase composition ofN: at 298.15 K 
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5.3.1 Effect of feed position 

Figure 5.41 shows the effect of feed position on the product purity. In the study, the total 
length of the bed has been divided into 209 equispaced grids The optimum purity is 
obtained on introducing feed at 60'*’ grid (60 / 209 grid) from extract recycle end of the 
bed. Figure 5.42 shows the final gas phase concentration profile in the bed undergoing 
feed step from different positions along the length of the bed. 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Dimensionless feed position from extract rcycle end 

Fig.5.41 Variation in product purity with different position of feed; F=60 mmol/half cycle, 

Ph= 4 atm, Pl=1.5 atm, Rr=1.0 



'ig.5.42 Final gas phase concentration profile ofN 2 in bed undergoing feed; F=60 mmol/half cycle, 

Ph=4 atm, Pl= 1.5 atm, Rr=1.0 





Here onwards we have taken 60"’ grid (optimum) from extract recycle end as the feed 
position (with total number of grids equal to 209). 

5.3.2 Original duplex PSA 

Figure 5.43 shows the mass balance sheet of the Original duplex for 60 mmol/half cycle 
of feed. The purity of extract (rich in N 2 ) and that of raffinate (rich in O 2 ) is only 72%. 
Figures 5.44 (a, b, c, d, e, f, g, h) show the gas and the solid phase concentration profiles 
in the bed undergoing feed, blowdown, purge and pressurization step respectively. 

It can be seen that, the gas and the solid phase concentration profiles in the bed 
undergoing feed is not stretched like that of CO 2 -N 2 and CH4-H2 systems. Subsequently, 
this bed will go for other steps like blowdown, purge and pressurization affecting the 
overall performance of the Original duplex. 

It has been observed that in MV-II duplex PSA, the bed profiles gets stretched throughout 
the length during the feed step leading to dramatic improvement in product purity. 




Unit of F, Et R, M and n - mmol 
Y - mole fraction of N 2 








































Pressurization 



Fig.5.44 (g) Gas phase concentration profile in the bed undergoing pressurization step (Original duplex) 
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Fig.5.44 (h) Solid phase concentration profile in the bed undergoing pressurization step (Original duplex) 
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5.3.3 MV-II duplex PSA 

P'igure 5.45 shows the mass balance sheet for MV-II duplex PSA. The purity of extract 
(rich in N 2 ) improved from 72% to 99.3% with 99.7% recovery and that of raffinate (rich 
in O 2 ) from 72% to 99.6% with 99.0% recovery over the Original duplex. 

The total gas coming out from the bed undergoing intermediate-blowdown step is very 
large (123 mmol) compared to that of CO 2 -N 2 and CH 4 -H 2 systems (<10 mmol). This is 
because of the light component, O 2 is also desorbed significantly when pressure is 
reduced in the bed from Ph to Pl unlike N 2 (non-adsorbing) in the CO 2 -N 2 and H 2 (non- 
adsorbing) in the CH 4 -H 2 systems. 

As a large amount has to be drawn from the bed undergoing blowdown to reduce the 
pressure, energy requirement is also large ( 12.2 kJ/mol of feed) compared to that of CO 2 - 
N 2 system (5.4 kJ/mol of feed) and CH 4 -H 2 (4.6 kJ/mol of feed). 

Figures 5.46 (a, b, c, d, e, f, g, h, i, j) show the change in the gas phase and the solid 
phase concentration with time in the bed undergoing feed, intermediate-blowdown, final- 
blowdown, purge and pressurization steps for the mass balance sheet given in Figure 
5.45. 
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Feed: ^igures 5.46 (a) and 5.46 (b) show the gas and the solid phase concentration 
profile of the feed step. Feed is introduced at 0.28 dimensionless bed length from the 
extract recycle end (top end). The Na from extract recycle and feed gets adsorbed and the 
ralfinate product, rich in Oa is drawn from the bottom. At the feed position, sudden 
change in the bed profiles can be observed as feed mixes with a stream of different 
composition. The bed profiles are stretched throughout the length unlike that of the 
Original duplex. 



Fig.5.46 (a) Gas phase concentration profile in the bed undergoing feed step (MV-II duplex) 



Fig.5.46 (b) Solid phase concentration prpfile in the bed undergoing feed step (MV-II duplex) 
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Intermediate-blowdown: Figures 5.46 (c) and 5.46 (d) show the gas and the solid phase 
concentration profiles of the intermediate-blowdown step. The solid phase concentration 
profile shows that as the pressure reduces from Ph (4 atm) to Pi (1.5 atm), the N 2 gets 
desorbed and then adsorbed like CO2 in the CO2-N2 system and CH4 in the CH4-H2 
system. The overall effect is that a little amount of N 2 (0.44 mmol) comes out from the 
bed. It can be concluded from the three systems studied that the locking effect is 
independent of whether only one component is adsorbed or both are adsorbed. The 
component locked inside the bed is always heavy component. 



Fig.5.46 (c) Gas phase concentration profile in the bed undergoing intermediate-blowdown step 


(MV-II duplex) 



Fig.5.46 (d) Solid phase concentration profile in the bed undergoing intermediate-blowdown step 

(MV-II duplex) 
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Final-hlowdown: I-igures 5.46 (e) and 5.46 (f) show the gas and the solid phase 
concentration proilic.s of the fmal-blowdown step. On evacuating the bed from top, the 
outgoing gas (extract product) is enriched with 99.29% N 2 .To draw 30 mmol extract 
product, the bed has to be evacuated from Pi (1.5 atm) to Pl (1.29 atm). 



N 2 enriched 
extract 


0.3 0.4 0.5 0.6 0.7 
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Fig.5.46 (e) Gas phase concentration profile in the bed undergoing final-blowdown step (MV-II duplex) 
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Fig.5.46 (f) Solid phase concentration profile in tbe bed undergoing final-blowdown step (MV-II duplex) 
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Purge: Pigures 5.46 (g) and 5.46 (h) show the gas and the solid phase concentration 
profiles of the purge step. The bed undergoing purge is at final-desorption pressure (1.29 
atm). As the raffinate recycle stream rich in O 2 enters the bed, the gas phase composition 
of N 2 decreases. This results in the decrease in the partial pressure of N 2 as the total 
pressure is constant. So, the N 2 gets desprbed and purged off from the bed. 
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Fig.5.46 (g) Gas phase concentration profile in the bed undergoing purge step (MV-II duplex) 
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Fig.5.46 (h) Solid phase concentration profile in the bed undergoing purge step (MV-II duplex) 
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Pressurization: Figures 5.46 (i) and 5.46 (j) show the gas and the solid phase 
concentration profiles of the pressurization step. The bed is pressurized from Pi (1.29 
atm) to Ph (4 atm). The composition of N 2 in the gas phase decreases as the inlet stream 
is rich in 99.6% O 2 . The amount required to pressurize the bed is very large (140 mmol) 
compared to that of the CO 2 -N 2 and CH4-H2 systems (< 10 mmol) because O 2 entering 
the bed is also getting adsorbed significantly unlike N 2 and H 2 . 



Fig.5.46 (i) Gas phase concentration profile in the bed undergoing pressurization step (MV-II duplex) 
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Fig.5.46 (j) Solid phase concentration profile in the bed undergoing pressurization step (MV-II duplex) 
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5.3.4 Compari.son of MV-II duplex and Original duplex PSA 

Figure 5.47 shows the comparative performance of the MV-II duplex and Original duplex 
PSA with variation in feed rate. It can be seen that the product purity in MV-II duplex 
improved drastically over the Original duplex in the feed range studied (60 mmol/half 
cycle to 1 00 mmol/half cycle). 

Effect of feed, F 



20 40 60 80 100 120 


Feed (mmol/half cycle) 

Fig.5.47 Effect of feed; PH=4atm, Pi =1.5 atm (MV-II), Pl =1.5 atm (Original) Rr=0.5 
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5.4 Duplex TSA 

Carbon dioxide-nitrogen system has been taken to study the possibility of separation in 
the duplex TSA. The parameters used in simulation are listed in Table 5.4. 


I'ablc .5 , 4 Parameters used in simulation of fractionation of carbon dioxide-nitrogen mixture (TSA) 


1 . 

2 . 

3 . 

4. 

5. 

6 . 

7. 

8 . 

9. 

10 . 

11. 

12 . 

13. 

14. 


Molar feed composition 

20% CO 2 , 80% N 2 

Feed position 

middle of the bed 

Adsorbent 

Zeolite 13X 

Particle diameter, cm 

0.0707 

Bed voidage 

0.4 

Bulk density, kg/m^ 

750 

Bed length, cm 

100 

Bed diameter, cm 

2.5 

Langmuir constant for CO 2 (bo), m /mol 

303.9 

Saturation constant for CO 2 (qso), mol/m^ 

22352 

LDF constant for CO 2 , s'^ 

0.808 

Heat of adsorption, J/moI 

Temperature parameters, IC' 

3.5E04 

k, 

-0.02816 

k2 

-0.02352 

Cycle time, s 

40 

Cooling (tc), J 

20 

Heating (tn), ^ 

20 


The temperature dependent Langmuir isotherm data for carbon dioxide-nitrogen mixture 
over Zeolite 13X reported by Cho et al. (1995) has been used. Nitrogen is non-adsorbing 
component. Figure 5.48 shows the adsorption isotherm of carbon dioxide on Zeolite 13X 
at different temperatures and pressures. 
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Fig.5.48 Adsorption isotherm of carbon dioxide on Zeolite 13X 

In duplex TSA, one bed undergoes cooling (adsorption) and the other bed undergoes 
heating (desorption). In the present study we have simulated the cooling step and heating 
step separately. The gas and the solid are at the same temperature. The initial gas phase 
concentration profile and temperature profile in the bed has been taken to be linear. 
During the steps, the temperature shift between the two ends of the beds has also been 
taken as linear. This is the ideal case, and has been taken to study the possible mechanism 
of separation. 

Cooling step 

Figure 5.49 shows the mass balance in cooling step. There are two inlet streams. Feed (10 
mmol) is introduced in the middle of the bed. The inlet stream from top of the bed has 
been taken as pure CO 2 (20 mmol). Figures 5.50 (a, b, c) show the temperature, gas phase 
concentration and solid phase concentration profile in the bed. The CO 2 entering from top 
end gets adsorbed in the first half. It can be seen from the concentration profiles that at 
feed position CO 2 gets desorbed. This may be because of the decrease in gas phase partial 
pressure of CO 2 as initial composition of CO 2 in middle of the bed is 0.5 and feed 
composition is only 0.2. 
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Temperature (K) 










Heating step 

Figure 5.51 shows the mass balance in heating step. The inlet stream from bottom of the 
bed has been taken as pure N. (20 mmol). Figures 5.52 (a, b, c) shows the temperature, 
gas phase concentration and solid phase concentration profile in the bed. The CO 2 gets 
desorbed as temperature inside the bed increases and then purged off from the top end. 
The composition of CO 2 in the outgoing stream is 71.8%. 

n = 68.3 
Y= 0.718 


MA,r 697.2 
MA.f= 648.1 
Msr 3.1 
Mg.p 3.9 
Ujr 700.3 
Mir 652.0 


n = 20.0 
Y= 0.000 

Fig. 5.5 1 Heating step in duplex TSA 
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Fig. 5.52 (a) Temperature profile in the bed undergoing heating step 
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Dimensionless bed length 

Fig. 5.52 (b) Gas phase concentration profile in the bed undergoing heating step 



Dimensionless bed length 


Fig. 5.52 (c) Solid phase concentration profile in the bed undergoing heating step 


The above results show that there is a possibility to get separation with the proposed 
duplex TSA cycle. However, further studies required to develop a feasible process. 
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Chapter 6 


Conclusions 

The Original duplex was not able to give sharp separation. To overcome this, 
modifications in the duplex PSA has been proposed to improve the performance. The 
parametric studies such as the effect of feed inlet position, raffinate recycle ratio, gas 
flow rate and desorption pressure were performed to assess the comparative performance 
of the Original duplex, MV-I duplex and MV-II duplex PSA with CO 2 -N 2 system (N 2 
was non-adsorbing), CH 4 -H 2 system (H 2 was non-adsorbing) and N 2 -O 2 system (both the 
components were adsorbing). 

In the original duplex PSA, the maximum purity obtained of the heavy component (CO 2 
and CH 4 ) was 90% and that of the light component (N 2 and H 2 ) was 95%. In the N 2 -O 2 
system, maximum purity of both the products obtained was only 75%. The MV-I duplex 
showed a marginal improvement in the product purity over that of the Original duplex 
whereas in the MV-II duplex, the purity of both the products improved dramatically to 
99+% with 99+% recoveries. The productivity in MV-II duplex PSA was 5-10 times 
greater than that of the conventional PSAs because of the small cycle time and extended 
mass transfer zone in the bed. 

The equalization and pressurization from both end of the bed viz. rich end (mode-I) and 
lean end (mode-II) was studied with the CO 2 -N 2 system. It was observed that the product 
purity for both the modes was about the same but the energy requirement in the mode-I 
was 1 .5 times more than that in the mode-II. 

The energy requirement with the N 2 -O 2 system (12.2 kJ/mol of feed) was found to be 
quite large compeired to that of the CO2-N2 system (5.4 kJ/mol of feed) and the CH4-H2 
system (4.6 kJ/mol of feed). This was because of the light component; O 2 was also 
desorbed significantly when pressure in the bed was reduced from Ph to Pi unlike N 2 in 
the CO2-N2 and H2 in the CH4-H2 system. 

Exploratory studies have been made to examine the feasibility of temperature swing 
adsorption. However, further studies required to establish the feasibility. 
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6.1 Scope of future studies 

The possible areas of future studies are as follows: 

> Fractionation of multicomponent mixture 

> Development of a duplex TSA cycle 

> Experimental validation of the simulation study 
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Appendix - A 


The governing partial differential equations (PDE’s) with the boundary conditions 
derived in Chapter-4 are discretized in space domain using central difference scheme of 
the finite difference method (FDM). 

A schematic diagram of the bed with grids is shown below. The total length of the bed 
has been divided into n+1 grids. Grid ‘0’ and ‘n+2’ represents ghost points. In the 
discretized equations at grid ‘1’ and ‘n+l’, we need values of variables at ghost point ‘0’ 
and ‘n+2’ respectively which are not known. These ghost point values are replaced by the 
known values at the nearby grids using the appropriate boundary conditions. 
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i i ! Point 
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The discretized governing equations and the boundary conditions for the duplex PSA and 
the duplex TSA are given below. In equations, subscript i represents grid and j represents 
component. 


A.l Discretized governing equations 
A.1.1 Duplex PSA 

Pressure profile; Equation (4.16) 
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Concentration profile in the gas phase; Equation (4.17) 
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Concentration profile in the solid phase; Equation (4.18) 
dq. 


dt 




A.1.2 Duplex TSA 

Pressure profile; Equation (4.14) 
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Concentration profile in the gas phase; Equation (4.15) 
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Concentration profile in the solid phase; Equation (4.18) 

dq,j 


dt 
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Temperature profile in the gas phase; Equation (4.21) 
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Temperature profile in the solid phase; Equation (4.23) 
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A.2 Discretized boundary conditions 
A.2.1 Duplex PSA 
> Blowdown step 
• at z = 0 
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> Purge step 
• at z = 0 


dz 

dz 
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= u 


yn+2 ~ yn 
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> Pressurization step 
• at z = 0 


dz 
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> Feed step 
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• at z = L 
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z; and z; represents the bed position just before and after the feed inlet position 
respectively. Expression for composition of the mixed stream, y,„ is given below. 


127 



From component mole balance at the mixing point we get, 
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where is calculated from total mole balance at mixing point 
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In expression of is calculated by discretizing the Blake-kozeny equation 

using backward difference scheme as 
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A.2.2 Duplex TSA 
> Heating step 
• at z = 0 
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> Cooling step 
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Zy and z/, represents tlie bed position just before and after the feed inlet position 
respectively. Expression lor composition and temperature of the mixed stream, and 
T^ ,„ is given below. 

From component mole balance at the mixing point we get, 
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where v,^ „ is calculated from total mole balance at the mixing point as. 
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In expression of v^- is calculated by discretizing the Blake-kozeny equation 

using backward difference scheme as 
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Appendix - B 


Integration scheme 

For numerical integration, Newton-Cotes fonnula (degree of polynomial - 4) has been 
used which is given as (Joe D. Hoffman, McGraw-Hill, 1 992) 

/ =-^A(7X+32^ + 127 ; +32/ +14/ +32/ + 12/ +32/ +14/ + 

+32/„+12/,+32/ + 7/.,) 

where n+\ is the number of grid points and h is the length of the equally spaced grid 
element. Total number of grid elements, n must be integral multiple of four. 
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